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. Abstract 

pi i' We explore the feasibility of the discrete flavor symmetry to explain the 

I pattern of neutrino masses and mixings. In the flavor basis, the neutrino mass 



matrix is taken to be invariant under ^3 symmetry at the zeroth order. The 
effects of breaking 5*3 symmetry by matrices invariant under different 52 sub- 
groups of S3 are studied. The resulting retrocirculant perturbation matrix leads 



> 

cn 

00 ■ to the perturbed 5*3 neutrino mass matrix having a trimaximal eigenvector as 

m ■ 

I suggested by the solar neutrino data. It is found that after the third order 

I perturbation the neutrino mixings only depend on the perturbation parameter 

and are consistent with the current experimental data leading to very small 
deviations from tribimaximal mixing. These perturbations in the S3 invariant 
^ ' neutrino mass matrix result in interesting interplay between the solar and the 

I reactor neutrino mixing angles. We also get a strongly suppressed range of 

effective Majorana mass which lies well within the reach of the ongoing exper- 
iments. 



1 Introduction 



The recent neutrino oscillation data have greatly enriched our knowledge of the flavor 
structure of leptons. After the discovery of neutrino oscillations, there has been a 
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considerable progress in determining the values of the neutrino mass-squared differ- 
ences and mixing angles relating the mass eigenstates to the flavor eigenstates. The 
remaining challenges are to determine the absolute mass scale of neutrinos, the mass 
hierarchy in the neutrino sector (sign of Am2^) and the value of 1-3 mixing angle (in 
particular if it is zero or not) amongst others. One has to understand the systemat- 
ics of the experimental data to unravel symmetries of the lepton mass matrices. A 
successful phenomenological Ansatz for the neutrino mixing matrix consistent with 
the present neutrino data was proposed by Harrison, Perkins and Scott [T] known as 
tribimaximal mixing (TBM) and is given by 



TBM 



( ^ ^ 

1 1 -1 



\ 



(1) 



\ ^ ^ ^ / 

\ V6 V3 V2 / 

Recently, non-Abelian flavor symmetries [21 [Sj [5] have been invoked intensively to 
account for TBM. These non-Abelian flavor symmetries can explain various, appar- 
ently unrelated, phenomena in flavor physics thus providing hints for the new physics 
beyond the Standard Model (SM). The basic features of the neutrino masses and 
mixings can be traced back to a minimum realization of the smallest discrete non- 
Abelian group 5*3 which is the permutation group of three objects. A large number of 
papers [31 13] have presented detailed models based on 5*3 symmetry. The permutation 
matrices in the three dimensional reducible representation are 



5(1) 



^ 1 \ 
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U ij 
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^(123) 
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^(132) 
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(3) 



^oio\ /ooi\ /ioo\ 

^(12) =10 0, 5(1=^) =010, ^(''^ =001, (4) 

^001/ \iooy \oiOy 

where the matrices in each equation belong to the same class of 6*3. The TBM 
mixing matrix Utbm is one of the transformations which relates the given basis to 
the irreducible one thus implying a close connection between TBM and S3 symmetry. 
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The neutrino mass matrix Mi, is invariant under which requires 



0, 



(5) 



where 5* is any of the six permutation matrices given in Eqns. (2, 3, 4). In the present 
work, we consider the possibihty in which the neutrino mass matrix is dominated by 
a term with symmetry. We then consider possible perturbations that violate this 
symmetry so that the perturbed neutrino mass matrix {M^) still satisfies the 
invariant constraint [7] 



M', 



M' 



M' -M' ,(/3^7^5). 



(6) 



Such connections have been considered in the literature from phenomenological point 
of view. In this analysis, we take the charged lepton mass matrix to be diagonal. This 
simple choice of basis can be made as the physical mixing matrix does not depend 
on the choice of basis. The property of basis independence allows one to rotate the 
charged lepton mass matrix to the fiavor basis where it is diagonal. If a horizontal 
symmetry exists it must be simultaneously a symmetry of the charged leptons as well 
as the neutrinos before the gauge symmetry breaking. After symmetry breaking when 
the fermions acquire mass, the charged lepton and the neutrino mass matrices should 
be constrained by different representations of the symmetry group in order to have 
non-zero mixing. The full group before symmetry breaking is at least a product group 
of the symmetry groups of the charged leptons and the neutrinos. Here we consider 
S's to be the residual group in the neutrino sector. For the charged lepton sector, 
symmetry can be taken as the residual symmetry which yields non degenerate 
diagonal charged lepton mass matrix 



2 5*3 invariant neutrino mass matrix 

The neutrino mass matrix implied by the commutator in Eqn. (5) is given by 



M^ = al + bD. 



where 
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Figure 1: Angle a between two equal length vectors a and a + 36. 



where a and h are in general complex and D is called the democratic matrix. The 
neutrino mass matrix can be diagonalized by the TBM matrix Utbm- The eigenvalues 
of are a, a+3b and a respectively. Here, the two degenerate eigenvalues correspond 
to the mass eigenvalues mi and m^, contrary to the experimental data where mi and 
m2 have smaller mass difference. This problem was overcome by Jora et al. [5] by 
introducing a Majorana-type CP violating phase a between the complex vectors a 
and a + 36 which does not affect the usual neutrino oscillations but affects the rate 
of neutrinoless double beta decay. The phase a can be adjusted to ensure equal 
magnitude (but different directions) of a and a + 36, thereby, leading to degenerate 
mass spectrum for 5*3 invariant My at the zeroth order. In the present work, the 
complex plane is oriented in a manner so that the parameter 6 is completely imaginary. 
The complex vector a lies in the fourth quadrant and is given by 



and the allowed range of a is < a < vr. The complex parameters a and 6 can be 
written in terms of the real free parameter X as 



For the Majorana-type CP violating phase a = 180°, the parameter x becomes zero 
and this value of a is, hence, not allowed. Democratic charged lepton mass matrices of 
the form bD and diagonal neutrino mass matrices of the form al have been discussed 



a = a e 



ia/2 



(9) 



\a\ = X sec{a/2), 
2 

\b\ = -xtan(a/2). 



(10) 
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earlier [6] in the context of S3 symmetry. The possibihty a = has, also been 
discussed in the literature |7j . S3 flavor symmetry in the context of Type (I+II) see 
saw has been discussed in ref . [8] . A generic feature of this approach is a democratic 
charged lepton mass matrix and a diagonal neutrino mass matrix at the zeroth order. 
In contrast, we consider the most general S3 invariant neutrino mass matrices at the 
zeroth order in the flavor basis where the charged lepton mass matrix is diagonal. 

3 Perturbations from 5*3 

In this section, we study the effects of perturbations from S3 invariant neutrino mass 
matrix Mj,, by the elements of 5*3 group itself. We follow a simple and highly predictive 
approach to obtain the results of perturbations order by order such that the neutrino 
mass matrix Mj, is invariant under S3 at the zeroth order. Parametrizing higher 
order perturbations of decreasing importance in terms of a small real dimensionless 
perturbation parameter A, the perturbed neutrino mass matrix M'^ is given by 

Ml = M, + fi{X[SP] + X^[S'^] + X^[S'']), (11) 

where is invariant under S3 and S^, S'^ and 5*^' can be any of S^^'^\ S^^'^^ and 
S^'^^'^ such that p ^ q ^ r and fx is the real parameter with dimensions of mass 
and magnitude of the order of one in the units of the absolute mass scale. The 
perturbation matrices are given in Eqn. (4) and belong to the same class of 5*3. 
Each perturbation matrix given above is invariant under a different 5*2 subgroup of 
5*3. The advantage of this approach of perturbation is that there are only four free 
parameters. Another appealing aspect of this approach is that the S3 symmetry is 
broken by its elements in the same representation. We see that even after breaking S3 
by above matrices, the S3 invariant constraint (Eqn. 6) remains intact which leads to a 
trimaximal eigenvector for the perturbed neutrino mass matrix M'^. We follow the line 
of argument given in Ref. [5] to decide as to which of the three perturbation matrices 
5'(i2)^ 5'{i3) g^j^^ 5'*^^^'' is the favoured one to be the first order perturbation. Since zeroth 
order neutrino mass matrix has two degenerate eigenvalues, we use quantum 
mechanical degenerate perturbation theory. Due to the double degeneracy between 
the eigenvectors corresponding to mi and 1713, the tribimaximal mixing matrix Utbm 
given in Eqn. (1) is not the unique one to diagonalize My. Rather it is a more general 
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matrix UtbmR{4>i3) where R{4>i3) is the rotation matrix in the 1-3 plane given by 



COS 013 — sin 
1 
\ sin 013 cos 013 J 



(12) 



In the ordinary quantum mechanical perturbation theory, the first order correction to 
the i*'^ eigenvector is given by ^S^. However, in case of degenerate perturbation 

£/j H/j 

theory, when Ei = Ej, the numerator has to vanish which implies the following 
condition for the (1,3) matrix element: 



(i?(013)^^^SM'5^^™i?(013))l3 = 0. 



(13) 



Amongst the three perturbation matrices only S^"^^^ satisfies the above condition thus 
leading to Utbm for = 0. The solution = tt, on the other hand, just changes 
the sign of first and third column of tribimaximal mixing matrix Utbm- However, 
the other two solutions = f and = ^ do not lead to the desired TBM form. 
The other two matrices 5'^^^^ and S'^^^^ also do not lead to the desired form of TBM 
mixing matrix, thus, making 5"*^^^^ as the favoured perturbation at the first order. 
The neutrino mass matrix after the first order perturbation with S^"^^^ becomes 



^ a + b + ijlX 



V 



b b 

a + b b + /iX 
b + /iX a + b J 



(14) 



The eigenvalues of this mass matrix are 



mi — a + IjlX, 

1712 — a + 3b + A*A, (15) 
777.3 — a — /iX. 

Here, we see that |mi| = |m2| and Im^l < \mi\ when A is positive i.e. after the 
first order perturbation \mi\ and |m2| remain degenerate, while |m3| < |mi| thus 
implying that our perturbation is in the right direction as experimentally the mass- 
squared difference Am^2 much smaller than Amfg. The neutrino mixing matrix 
after this first order perturbation is still Utbm- The next step is to add the higher 
order perturbations to remove the degeneracy in the magnitudes of 777i, 7772 and also 
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to deviate from TBM. This can be done by taking X^S^^"^^ as the second order and 
A^S*^^^-* as the third order perturbation or X^S^^^^ as the second order and X^S^^'^^ as 
the third order perturbation. We investigated both these possibihties in our analysis 
and found that the choice of the second and the third order perturbation matrices 
only affects the predictions for the 2-3 mixing angle which we discuss later. Taking 
^2g{i2) ^Yi^Q second order perturbation and A'^S''-^^-' as the third order perturbation, 
the perturbation matrix becomes 



f A 

A2 
VA3 



A2 
A3 



A A2 y 



(16) 



Hence, the perturbed neutrino mass matrix becomes 



M' 



a + b + fiX b + jj,X^ 
b + jjX^ a + b + fiX^ 
\ b + jjX^ b + jjX 



b + jiX^ 
b + 12X 
a + b + /iA^ J 



(17) 



It is interesting to note that the perturbation matrix in Eqn. (16) is retrocirculant 
having a trimaximal eigenvector (■^, :^)^ O HO]. Also the perturbed neutrino 
mass matrix M'^ is a magic mass matrix (TUl [H] . The eigenvalues of the perturbed 
neutrino mass matrix upto the third order perturbation are given by 



mi ^ xsec(|)e + /i(A 



A3' 



m2 



xsec(f)e-^"/2 



2ixtan(f) +/i(A + A2 + A^ 



^ xsec(f )e-^-/2 + /i(-A + f + f). 



The magnitudes of the eigenvalues which are the neutrino masses are given by 



\mi\ 
\m2\ 
\ms\ 



xsecf [1 + cos2f{^ + 



-,2 ajXfM 

X 

xsecf [1 + cos2f{^ + 



xsecf[l + cos2f{f + 



A^(^^— x/^— /i^ cos^ ^) 

\'^(^l'^+2x^l-^? cos^ f ) 



\^(x'^ ^i+Ax^'^+i^i^ cos^ ^—Axfi^ cos^ #—4/^3 cos' 



r)- 



2x2 

x'^{^?+x^JL-^? cos^ |) 
2^ 



8x3 }] 

A3 (2x2^+2x^2-/^3 C0s2 f -2x/j2 cos2 |+/x3 COS* f ) 



2x3 
"^3" 



A3(x2^-4x/i2+4/i3 cos2 |+4x/t2 cos2 ^-if? cos* f ) 



}], 
}]• 



(19) 



These neutrino masses are made real positive by the phase matrix 



P 













(20) 
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where 



r 



a 



■tan(f) 



8 I 



p 



tan(f) 

2^"" Vi+M(A+A2+A3) 

-tan(|) 
l+^(_A+f+f) 



I tan ^( 



(21) 



)■ 



It can be seen from Eqn. (19) that the neutrino masses depend on the Majorana-type 
CP violating phase a. However, the lepton number conserving neutrino oscillations 
do not depend on this phase. Therefore, we find the mass-squared differences are 
independent of the Majorana phase. 



^2 



(22) 



Notice that the leading order term in the solar mass-squared difference Am^a is of 
the order of which makes it naturally small to be consistent with the experimental 
data. The neutrino mixing matrix after the third order perturbation is given by 



U' 



^33 I 



(23) 



where 




t^33 



(24) 
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An interesting feature of this type of breaking of Ss is that the neutrino mixing matrix 
is completely determined from the retrocirculant perturbation matrix and is, thus, 
independent of the parameters appearing in the zeroth order neutrino mass matrix 
Mj,. As a consequence, the deviations from TBM are entirely determined in terms 
of single unknown perturbation parameter A. The neutrino mixing angles upto the 
third order perturbation are given by 

1 X3 

smelo^- + — -—, (25) 
12 3 24 24' ^ ' 

sin^L^y-y, (26) 

and 

sinel,^ - + ---. (27) 

23 2 4 8 ^ ' 

The solar and reactor neutrino mixing angles are related by 

sin^^3 = 3sin0^^ - 1. (28) 

As pointed out earlier in this section, the second and third order perturbation matrices 
can be interchanged. As a consequence of this interchange, the predictions for all the 
parameters remain intact except for the 2-3 mixing angle 623 which is given by 

1 A 

sin^^^ ('29) 

23 2 4 8 ^ ^ 

It can be seen that the 2-3 mixing angle in this case is shifted to below maximal by 
the same amount as it was above maximal in the earlier case. 



4 Numerical results 

The above results are obtained when the perturbation parameter A is positive. Since A 
is real it can also be negative. When A is positive, we get the inverted hierarchical mass 
spectrum for neutrinos and a negative A gives a normal hierarchical mass spectrum. 
We use the constraints for the mass squared differences from the known oscillation 
data [12] viz. 

A™2 _ 7^7(+0.16,+0.34,+0.52) .^_5 -5^2 
^'''-12 — '■"'(-0.19,-0.36 -0.53) ^ 

Aml^ = ±2.39|^og lg 2o,lo.33)^ ^ 10 ^eV"^. 

(30) 



9 



For normal hierarchy i.e. for negative A, the parameters /i and A are constrained to 
he in the range (0.2 - 0.3)eV and ((-0.037) - (-0.095)) (Fig. 2(b)) respectively at 
the 3(7 CL. in this model whereas the free parameter x and the Majorana phase a 
are constrained to the range (0.035 — 0.09)eV^ and (0° — 160°) respectively at the same 
confidence level. The model, also, leads to small deviations of the mixing angles from 
TBM, given by 

^12 = 35.26° - 35.29°, (31) 
= 0.8° - 2.2°, (32) 
023 = 43.45° - 44.45°. (33) 

while the allowed ranges for these mixing angles [121 given by 

n _ qo Q(-o(+1.16,+2.43,+3.80) 
Pl2 — 'J'J-yD(_i.i2 -2.13 -3.10) ' 
a _ /IQ nc°(+4-18,+7.83,+10.32) 

013 < 12.38°(3a). (34) 



For inverted hierarchy i.e. positive A the allowed ranges for the parameters /i and 
A are (1.44 - 0.2)eV and (0.017 - 0.043) (Fig. 3(b)) respectively at the 3a C.L. 
whereas the free parameter x and the Majorana phase a are constrained to the range 
(0.2— 1.44)e\^ and (0°— 160°) respectively at the same confidence level. The deviations 
of the mixing angles are given by 

= 35.265° - 35.269°, (35) 
^13 = 0.35° - 0.90°, (36) 
023 = 45.25° - 45.61°. (37) 

It can be seen from the neutrino mixing matrix given in Eqn. (23) that there is no 
Dirac-type CP violation in this model as A is taken to be real. The CP violation 
in this model is only due to the phase a which gives Majorana-type CP violation. 
However, A may, in general, be complex and there will be another CP violating phase 
in the mixing matrix for complex A leading to Dirac-type CP violation. The possible 
measurement of the effective Majorana mass M^e in the neutrinoless double beta de- 
cay searches will provide information on the Majorana-type CP violating phase. The 
analysis of Mge is significant because a stringent constraint on the value of effective 
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Figure 2: Correlation plots for the normal hierarchy. 
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Figure 3: Correlation plots for the inverted hierarchy. 

Majorana neutrino mass Mee was obtained in the ^^Ge Heidelberg-Moscow experi- 
ment [13] |Mee| < 0.35eV. There are a large number of forthcoming projects such as 
SuperNEMO [14j, CUORE [15], COURICINO [16j and GERDA [T7] which aim to 
achieve a sensitivity below O.OleV. The effective Majorana mass which determines 
the rate of neutrinoless double beta decay is given by 



= \miU[l + maf/fs + m^U'^l \ . (38) 

The value of M^e for normal and inverted hierarchies are (0.015 — 0.083)eV^ and 
(0.04 — 0.12)ey. The correlation plots between the effective Majorana mass and the 
perturbation parameter A is given in Fig. 2(b) and 3(b). 
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5 Summary 



Wc studied the breaking of invariant neutrino mass matrix in the flavor basis by 
perturbing it to the third order of a small perturbation parameter to examine its 
consequences for neutrino masses and mixings. The neutrino mass matrix at the 
zeroth order gives two degenerate eigenvalues mi and thus leading to results 
contrary to the experimental data. This problem is flxed by the introduction of the 
Majorana phase which renders the magnitude of all the three eigenvalues equal at 
the zeroth order. This degeneracy is broken by the introduction of perturbation 
matrices chosen to be 5*3 group matrices invariant under three different 5*2 subgroups 
of Ss- The higher order perturbations are suppressed by increasing powers of the 
small perturbation parameter. The quantum mechanical degenerate perturbation 
theory uniquely chooses the first order perturbation matrix while the second and the 
third order perturbations are simultaneously added. The mass-squared differences are 
found in terms of the perturbation parameter A, the free parameter x and parameter /i 
which is of the order one. Interestingly, the mixing matrix and the mixing angles in the 
model are independent of the parameters of the S3 invariant neutrino mass matrix and 
are determined only in terms of perturbation parameter A. This scheme of S3 breaking 
results in very small deviations from the TBM when confronted with experimental 
inputs of the mass-squared differences. Both normal and inverted hierarchical mass 
spectra are allowed in this model. The cause of CP violation in this scenario is the 
Majorana-type CP violating phase. The allowed range of the effective Majorana mass 
obtained in the present study lies well within the reach of the forthcoming experiments 
providing a stringent test of this model. 
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